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, Abstract. We prove that a Poisson-Newton formula, in a broad sense, is associ- 

ated to each Dirichlet series with a meromorphic extension to the whole complex 
^\ , plane. These formulas simultaneously generalize the classical Poisson formula and 

' Newton formulas for Newton sums. Classical Poisson formulas in Fourier analysis, 

classical summation formulas as Euler-McLaurin or Abel-Plana formulas, explicit 
formulas in number theory and Selberg trace formulas in Riemannian geometry 
appear as special cases of our general Poisson-Newton formula. We also associate 
to finite order meromorphic functions general Poisson-Newton formulas that yield 
, many classical integral formulas. 
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1. Introduction 

All classical Poisson formulas for functions in Fourier analysis result from the gen- 
eral distributional Poisson formula 
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which is an identity of distributions identifying an infinite sum of exponentials, con- 
verging in the sense of distributions, and a purely atomic distribution. This distribu- 
tional formula is related to the simplest finite Dirichlet series 

f(s) = 1 - e~ Xs . 

It is interesting to observe that on the left hand side of (OQ) we have an exponential 



sum 
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where the sum runs over the zeros p n = ^n, n G Z of /, and on the right hand 
side of (CQ) we have a sum of atomic masses at the multiples of the fundamental 
frequency A. One can say that the frequencies associated to the zeros are resonant 
at the fundamental frequencies. Taking the Fourier transform we obtain the dual 
Poisson formula that is of the same form where we exchange zeros and fundamental 
frequencies. Thus the fundamental frequencies are also resonant at the zeros. 

The main purpose of this article is to show that this is general and to each mero- 
morphic Dirichlet series / we can associate a distributional Poisson formula 

(2) W(f) = J2 n p ePt = E^'^'W , 

P k 

where the first sum of exponentials runs over the divisor of /, i.e., zeros and poles 
p with multiplicities n p , and the second sum runs over non-zero sequences k = 
(k\, &2, ■ ■ ■) € N°° of non-negative integers, all of them zero but finitely many, and 
(A, k) = ^2Xjkj. The equality holds in M*,. Conversely, we prove that any such 
Poisson formula comes from a Dirichlet series. 

The distribution 

p 

is well defined in and is called the Newton-Cramer distribution of /. We name it 
after Newton because it appears as a distributional interpolation of the Newton sums 
to exponents t G R, since in the complex varia bk3 z = e s the zeros are the a = e p so 

a 

and for integer values t = m G Z we get (in case of convergence) the Newton sums 

W(f)(m) = S m = J2^ m - 

a 

There is a precise theorem behind this observation. We show that our Poisson- Newton 
formula for a finite Dirichlet series / with a single fundamental frequency is strictly 
equivalent to the classical Newton relations. This is the reason why we name also 
after Newton our general Poisson formulas. 

Writing p = 27 we see that the sum W(f) of the left hand side of (|2J) is the Fourier 
transform of the atomic Dirac distributions <5 7 and we can formally write 

7 k 



The variable z = e s or better z — e s is the proper variable when dealing with Dirichlet series. 
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The form of this formula, relating zeros to fundamental frequencies, strongly reminds 
other distributional formulas in other contexts. In number theory, more precisely in 
the theory of zeta and L-functions, the same type of identities do appear as "explicit 
formulas" associated to non-trivial zeros of the zeta and other L-functions. These ex- 
plicit formulas, when written in distributional form, reduce to a single distributional 
relation that identifies a sum of exponentials associated to the divisor of the zeta or 
L-function and an atomic distribution associated to the location of prime numbers. 
Usually the sum runs over non-trivial zeros, and the sum over trivial zeros appears 
hidden in other forms as a Weil functional, which is classically interpreted as corre- 
sponding to the "infinite prime"Q For that reason, Delsarte labeled this formula as 
"Poisson formula with rest" (see [TTJ), the "rest" refers to the sum over the trivial 
part of the divisor. More precisely, for the Riemann zeta function, we have in M.+ 

J2 n P ePt + w o(f) = J2J2 lo s pSki °zp ' 

p p k>l 

where the sum on the left runs over the non-trivial (i.e., non-real) zeros p, and the sum 
over p runs over prime numbers. Conjecturally, the non-trivial zeros are simple, i.e., 
n p = 1. The term Wo(f) is the sum over the trivial (real) divisor and is computable 

W (f)(t) = -e* + ^V*" = + -2TTT ' 

n>l e 

and corresponds to Delsarte "rest", or to the Weil functional of the infinite prime. 
Also we have in this case 

p 

where 

v(t) = , 

»7>0 

is the classical Cramer function, studied by H. Cramer [9], where p = \ + ij. This 
motivates that we name our distribution W(f) also after Cramer. 

In Riemannian geometry, we have the same structure for the Selberg trace formula 
for compact surfaces with constant negative curvature. With the relevant difference 
that Selberg zeta function is of order 2, which gives a "rest" of order 2 also. Selberg 
formula relates the length of primitive geodesies, which play the role of prime num- 
bers, and the eigenvalues of the Laplacian, which give the zeros of the Selberg zeta 
function. For non-negative constant curvature, the formulas are of a different nature 
and the distribution on the right side are no longer simple atomic Diracs, but also 
higher order derivatives appear. This will be discussed elsewhere. In the context of 



'It may be more appropriate to talk of the prime p = 1. 
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dynamical systems and semiclassical quantization, we have Gutzwiller trace formula, 
which relates the structure of the periodic orbits of a classical mechanical system to 
the energy levels of the associated quantum system. 

The interpretation and analogy of these formulas with "Poisson formulas" was 
noticed long time ago. We should mention in particular the classical work of A. P. 
Guinand [13], J. Delsarte [IT] . A. Weil [28J, and results related to Hamburger theorem 
[15] , [19] . Already the title of Delsarte's article points to the Poisson flavor of these 
formulas "Formules de Poisson avec reste". More recently this analogy between 
Poisson and explicit formulas and its relation with zeta-regularization is studied for 
the Selberg trace formula for surfaces with constant negative curvature by P. Cartier 
and A. Voros in [7|. General Poisson formulas for Riemannian manifolds relating the 
spectrum of positive elliptic operators and the length spectrum were developped by 
J. Chazarain [10] and J.J. Duistermaat and VW. Guillemin [12], and the Dirichlet 
series associated to the spectrum of the heat equation by S. Minakshisundaram and 
A. Pleijel [21], after foundational work by T. Carleman [8]. 

Our goal is to put in the proper context, generalize and make precise the analogy 
of Poisson and trace formulas, and derive a general class of Poisson formulas that 
contain all such instances. More precisely, to each meromorphic Dirichlet series of 
finite order we associate a Poisson-Newton formula. All relevant known formulas can 
be generated in this way. On the other hand the fact that explicit formulas in number 
theory and Selberg trace formula can be seen as a generalization of Newton formulas, 
seems to be a new interpretation. 

It is important to remark that in our general setting the Poisson-Newton formulas 
are independent from a possible functional equation for the Dirichlet series /, contrary 
to what happens in classical formulas. As a matter of fact, we do associate a Poisson- 
Newton formula to any Dirichlet series with no functional equation. This is sometimes 
hidden in the classical theory where explicit formulas and functional equations come 
hand to hand. For instance, most proofs in the literature derive the explicit formula 
for the zeta function using its functional equation, and the "rest" term, borrowing 
Delsarte terminology, is computed from the fudge factor from the functional equation. 
Nevertheless our approach shows that the functional equation is not related to the 
existence of an explicit formula. Moreover, for an / having a functional equation, the 
"rest" term in our Poisson-Newton formula emerges from the non-symmetric part of 
the divisor of /. 

Although independent, these questions are interrelated. As is well known, it is 
a classical and basic procedure since Riemann foundational memoir [23] to derive 
functional equations for the Dirichlet series of / from Poisson formulas for other 
Dirichlet series (for the function in the case of the Riemann zeta function). 
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We also derive a general Poisson-Newton formula associated to finite order mero- 
morphic functions f(s) which are not necessarily Dirichlet series, but have their divisor 
contained in a left half plane. In this general situation, the Newton- Cramer distribu- 
tion is still defined by the exponential series (converging in R + as distribution) 

p 

This time, the distribution is no longer a sum of purely atomic measures in R^. A 
particularly important case is when the divisor of / if left-oriented, i.e., contained 
in a left cone. Then the Newton-Cramer distribution is a ^-distribution and it is 
a distribution given by an analytic function in R^. An application of our general 
Poisson-Newton formula gives a collection of classical formulas: Gauss formula for 
the logarithmic derivative of the T-function, Binet formula for the logarithm of the 
T-function, general Gauss and Binet formulas for higher order Barnes T-functions, 
etc. 

The structure at of the distributions appearing in the Poisson-Newton formula is 
interesting. In the construction of Newton-Cramer distribution we have some param- 
eter freedom that is irrelevant for the structure of the distribution in R*, but not at 
0. But precisely the variation of this parameter gives a distributional form of differ- 
ent infinite Euler-MacLaurin type formulas associated to each Dirichlet series. More 
precisely, the classical Euler-MacLaurin formula (as well as Abel-Plana summation 
formula) can be derived from the simplest case of the Dirichlet series f(s) = 1 — e~ s . 
The distributional infinite Euler-MacLaurin formula sheds some light on Ramanujan's 
theory of the "constant" of a diverging series. Most of the other summation formulas 
over the semi-group generated by the frequencies of the Dirichlet series seem new. 

2. Dirichlet series 

We consider a non-constant Dirichlet series 
(3) /( S ) = l + ^a n e- A "% 

n>l 

with a„eC and 

< Ai < A 2 < • • • 

with (A n ) a finite set (equivalently, take the sequence (a n ) with all but finitely many 
elements being zero) or A n — > +oo, such that we have a half plane of absolute con- 
vergence (see [18] for background on Dirichlet series), i.e., for some a G R we have 

\a n \ e~ Xn ° < +oo . 

n>l 
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It is classical ([IB], P-8) that 

_ log(|ai| + |a 2 | + . . . + \a n \) 
a = lim sup . 

The Dirichlet series ([3]) is therefore absolutely and uniformly convergent on right 
half-planes > er, for any o > a. 

We assume that / has a meromorphic extension of finite order to all the complex 
plane s 6 C. We denote by (p) the set of zeros and poles of /, and the integer n p is 
the multiplicity of p (positive for zeros and negative for poles, with the convention 
n p = if p is neither a zero nor pole). The convergence exponent of / is the minimum 
integer d > such that 

^2 Kl \p\~ d < +°° • 

p+Q 

We have d = if and only if / is a rational function, which cannot be a Dirichlet 
series, thus d > 1. Indeed we always have d >2 (see Corollary 13.81) . The order o of 
/ satisfies d < [o] + 1. 

Since / has finite order, we have the Hadamard factorization of / (see [2], p. 208) 

f( 8 ) = s no e^Y[E m (8/p) n <' , 

where m = d — 1 > is minimal for the convergence of the product with 

E m (z) = (1 - z)e* + ** 2+ - + ^ m , 

and is a polynomial uniquely defined up to the addition of an integer multiple of 
2ni. The genus of / is defined as the integer 

g = min(degQ/,m) , 

and in general we have d < g+1 and g < o < g + 1 (see [2], p. 209). For a meromorphic 
Dirichlet series we prove that in fact d = g + 1 (see Corollary 13. 7\i . 

The origin plays no particular role, thus we may prefer to use Hadamard product 
with origin at some a G C, 



f(s) = {s- a)^e^^l[E m (j^)^ 



We have, uniformly on 3?s, 

lim f(s) = 1 , 

3ffs->+oo 

thus 

<ti = sup < +oo 
p 
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so f(s) has neither zeros nor poles on the half plane 9fts > o\. Sometimes in the 
applications o\ is a pole of / because when the coefficients (a n ) are real and positive 
then / contains a singularity at a by a classical theorem of Landau (see [18], Theorem 
10, p. 10). The singularity is necessarily a pole by our assumptions, and in general 
ox = a. 

Associated to the divisor div(/) = J2 n pPi we define a distribution W(f) = 
^n p e pt on W + . We do this as follows. 

Consider the space S of C°°-functions of rapid decay on R (i.e., <p G S if and only 
if for any n, m > 0, \t\ n D m ip — > 0, as t — > ±oo). The dual space S' is the Schwartz 
space of tempered distributions. As V = C S, we have that S' C V, where V is 
the space of distributions. 

Lemma 2.1. For finite sets A, consider the family of locally integrable functions 

\peA J 

There is a family of distributions Wa{J) which coincides with W / a(/) in M*, and which 
converges in K (over the filter of finite sets A), to a distribution W(f) in V . 

This distribution has support contained in R + , is Laplace transformable, and e~ ait W(f) G 
S' . More precisely, if o~\ is not a zero nor pole (resp. it is a zero or pole), e~ ait W(f) is 
the (distributional) d-th derivative of a uniformly bounded continuous function (resp. 
continuous function) on K with support in R + . 

More precisely, we have 

W(f) = e^^((K d (t)-K d (0))l R+ ) , 

where 

Proof. We prove first the lemma when <j\ is not a zero nor pole. We define 
(4) K l(t) = £ {j^-,^') V- 

Then for £ > d, Kg is absolutely convergent for t G 1R+ since 

| e O-<Ti)t| _ e 3t(p-ai)t ^ j 



8 V. MUNOZ AND R. PEREZ MARCO 

and Ki is a uniformly bounded function in R, continuous in R*, since 



< OO . 



p - <7i| 



The function 

F,(t) = (^(t)-^(0))1 M+ 
is a uniformly bounded continuous function on R, for t > d. 

For a finite set A, denote by 



\ ( — (T-i ; ; / 

peA 

the corresponding sum over p <E A, and -^,a(0 = (Ke,A(t) — -fQ,A(0))l]R + . On R*, 

= (5>^ ) ^+ = ^(*)- 

We consider 

W A(f) - e^ d F d , A , 

taking the distributional derivative. 

For a smooth function (resp. function with polynomial growth) K on R and a test 
function if with compact support (resp. in the Schwarz class) we have 

D -(Kl R+ ),v) = -(Kl R+ ,v') 



Dt 

r+oo 

K(t)<p'(t) dt 



r+oo 

= -[K(t)<p(t)]+°° + / K'(t)<p(t)dt 
Jo 

= K(p)<p(0) + (K'lv L+ ,<p) , 



thus 



^- t (Kl R+ ) = K'1 R+ + K(0)5 . 



Then since K' e A = K e _ ljA we get 
—.F^A = K 0tA (t) + K hA (0)5 + K 2 , A (0)8' + ... + K d _ M (0)#- 2) 



Dt d 

^ (i)+ fe^)* + (S^)* + ''' + (S^ 



(d-2) 
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Thus the difference between W A (f) and W A (f) is a distribution supported at {0}. 

We have the convergence F dA — > F d , uniformly as continuous functions on R. Thus 
we have the same limit F d>A — > F d in the distributional sense. Then taking the limit 
as distributions, W A (f) — > W(f), where 

W (f) = e^ d F d , 

which is the d-th derivative of a uniformly bounded continuous function on R with 
support on R + , as stated. 

When o"i is part of the divisor, then we do the same proof with 



•+' 



which adds to W(f) a term n ai e ait . □ 

Note that we can write 

W(f)\ K = ]imW A (f)\ K = J2 n P ePt > 

p 

as a distribution on Rl. But if d > 2, the family of distributions W A (f) is not 
converging to a distribution in R because the sums 

y n P 

are not absolutely convergent for £ — 1, . . . , d — 1 (by the definition of d). On the 
other hand, the same argument shows that for £ > d, Ke t A has a limit Kn in the sense 
of distributions, and for any k > 

D k D k 



Dt k ^ Dt k 

It is important to note that W(f)\u.* + is independent of the choices of <j\ and of 
taking d larger than the exponent of convergence. The only dependence on <j\ and d 
is located at the structure of the distribution at 0. 

Proposition 2.2. Let d be the exponent of convergence of f as before, and let d! > d. 

We define for a e C — {p}, £ > 0, and a finite subset A C {p} 



F ttA (t, a) = (K ttA (t, a) - K ttA (0, cr))l n 
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then the following limits exist in the sense of distributions 

W(f,d',a) = \jme*j^K d , tA (t,<r), 

W(f,d',a) = Hme^ t ^ F F d , iA (t,a), 

and W(f, d', a) |r* = W(f, d', a) |r* = W(f) |r* is independent of a and d! > d. 

Proof. We shall deal with the first case, the second one is similar. The existence of 
the limit is proved as before. Take 

W A (f,d',a) = e^^K^A^o-) = e*^ [e^-^K A , A {t^)) , 

where <j\ is the one considered before. Now, since d! > d, K d r jA (t,ai) converges to a 
distribution K d >(t, <7i). So W A (f,d',a) converges as distribution to 

W(f,d',a) = e*^ (e^-^KAt,^)) = e'^K^a). 

In M.*_ the independence on a is clear since the distributions vanish. In first 
we note that for £ > 0, 

d 
at 

and 

^(Kl m+ ) = K'l R+ + K(0)5 . 
So we have, in the sense of distributions: 
j^rKe,A = K e - 1)A + K i<A (0, a) S 

§g(Ki,A) = K , A + K hA (0, a) 5 + K 2fA (0, a)5' + ... + K tA (0, a) 5^ 

Using this, we get 

jjd> d ' 
(5) e at -^K d , A (t, a) = e^K 0A (t, a) + e CTi ^ K t , A (0, a) . 

i=i 

From this last expression we see that away from the distributional limit, that we 
know to exist, is independent of a and d', since e at K 0tA (t, a) = J2 peA e pt is indepen- 
dent of a, and the other summand is supported at zero. □ 

In section [4] we study in more detail the parameter dependence at 0. 
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3. Poisson-Newton formula 

On the half plane 9fts > o~x, log f(s) is well defined taking the principal branch of 
the logarithm. Then we can define the coefficients (6k) by 



(6) - log f{ S ) = - log ( 1 + a n Z~ XnS J = E 6k e 

V n>l / keA 



-<A,k)s 



where A = {k = {k n ) n >i \ k n £ N, ||k|| = ^ |/c n | < oo, ||k|| > 1}, and (A, k) = 
\iki + . . .+\iki, where k n = for n > I. Note that the coefficients (6k) are polynomials 
on the (a n ). More precisely, we have 

m ,, (-!)'"■" ||k||! yj kj 

Note that if the X n are Q-dependent then there are repetitions in the exponents of 
3.1. Hadamard interpolation. 

Lemma 3.1. Consider a discrete set {p} C C with the property that 

|n p | < +oo . 

Let (7i G C. VFe /iaue t/mt 

n„ , (s-ffi)''- 1 1 



S-Gi ' ^ P (p - (Ti)^ 1 S - p ' 

zs a meromorphic function in C, and /ias a simple pole with residue n p at each p. 

Proof. We start by noting that 

1 ^ (.s-a,)' !.s- ff| )"-' I 



(S — (Ti) 1 ' ' " • 

1=0 



- y 

p-s r^{p-(y\) l+l {p-a l ) dl p-s 



Consider a disk D(0, R) C C, and split G(s) = Gi(s)+G 2 (s), where Gi(s) corresponds 
to the sum of those p G D(0,R), and G^s) to the sum over the remaining p's. 
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Now, for s e D(0, R/2), we have 

(8) |G 2 (s)| < C\s - <7i| d 5] \n p \ \p\- d < oo , 

p 

thus we get the absolute and uniform convergence of the series in D(0, R/2). As Gi(s) 
has simple poles at p with residues n p , in D(0, R/2), we get the required properties 
for G(s) in D(0, R/2). This happens for every R > 0, thereby the result. 

□ 



Now we can define the Hadamard interpolation associated to the divisor ^2n p p. 
We define it up to a multiplicative constant which is irrelevant when we consider its 
logarithmic derivative as we will do. 

Definition 3.2. We define the Hadamard interpolation as 



f H (s) = exp (^J G(s) ds^j . 



The divisor of /# is Div f H = Yl n pP- 

Definition 3.3. Consider a meromorphic function f with divisor ^2n p p. Fix o\ as 
above. The discrepancy of f is defined as the difference of the logarithmic derivatives 

Ps -T«-J- G -J 

Lemma 3.4. The discrepancy Pf is a polynomial of degree < g — 1. 



Proof. We recall the Hadamard factorization 
f(s) = (s-a 1 ) n ^e Q f^ Yl 

where Q is a polynomial of degree < g which is uniquely defined up to the addition 
of an integer multiple of 27ri. The product is absolutely convergent because of the 
definition of the convergence exponent d. Taking the logarithmic derivative, we have 
that 

j = Q'f^ + G(s) , 

so Pf = —Q'j and the result follows. □ 



E<t-i ( '-^ ) 

\P-<7lJ 
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3.2. Poisson-Newton formula. The main result is the following Poisson-Newton 
formula associated to the Dirichlet series /. Consider its polynomial discrepancy Pf 
of degree < g — 1, 

P f (s) = c + cis + . . . + Cg^s 9 - 1 . 
The inverse Laplace transform of Pf is the distribution supported at {0} 

C-\Pf) = c 5 + Cl 5' + ... + c g - 1 8$'- 1) . 



Theorem 3.5. As distributions in R we have 

9-1 

^(/) = ^c fe # ) + ^(A,k)6 k 5 {XM) . 

k=0 keA 

The structure at of W(f) depends on the function / and its comparison with 
the Hadamard interpolation. The structure out of only depends on the divisor 
and is independent of parameter choices. In some sense it is the most canonical part. 
Sometimes we refer to the "full Poisson-Newton formula" the one of the main theorem 
with the structure at 0. 

Corollary 3.6. As distributions on R* + we have 

keA 

Proof. We prove the theorem by taking the right-sided Laplace transform of W(f) 
(we use the interval [— 1, oo)): 



roo M 

= jf (-l) d (K d (t) - K d {Q))—e^dt 

= n gl (-l)' , (<7l ~ 8)<< f + °°t d e^ t dt+ 
d- Jo 

n ai (s - o x 



+oo 

e (<Tl-« 



s — o\ P (p — a i) d \p — s °~i — s 



\d-l 



n °i _x^ n ( s - a iY 

s-oi p (p _ aiY' 1 p - s 

G(s). 
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On the other hand, consider the distribution 

k 

Its Laplace transform is 

(V, e"* s )[_ li00) = ( 5^(A, k) 6 k 5 (Ajk> , e" 



is 

/ 

[-1,00) 



k 



= -(-log/(a))' 

f(£) 
" /(«) 

= G( a )-P/( a ) 

= (W(f),e- st ) hhoo) -P f (s). 

By uniqueness of the Laplace transform for distributions (see [30], Theorem 8.3-1, 
p. 225), we have 

w(f) = v + c~\p f ), 

where C^ 1 {Pf) is the inverse Laplace transform of the polynomial Pf. This is a 
distribution supported at {0}. Hence we get the theorem and the corollary 

W{f)\ K = V. 

□ 



Just inspecting the order of the distributions appearing in both sides of the Poisson- 
Newton formula, we get an interesting corollary for Dirichlet series. We know that 
d < g + 1. In fact we do have equality d = g + 1. 

Corollary 3.7. For a meromorphic Dirichlet series we have 

d = g + l = degQ/ + l = degP/ + 2 . 

Proof. We inspect the order of the distributions in the Poisson-Newton-formula. We 
recall that W(f) is, as distribution, the d-th derivative of a continuous function. But 
8q^ is not the d-th derivative of a continuous function for I > d—1. Thus deg Pf < d—2 
so deg Qf < d — 1 and g < d—1, hence g + 1 = d. □ 

It is clear that d > 1 for a meromorphic Dirichlet series, but we have in fact d > 2. 
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Corollary 3.8. For a meromorphic Dirichlet series we have a convergence exponent 
at least 2 and order at least 1: 

d>2 

and 

o > 1 . 

Proof. As before we inspect the order of the distributions in the Poisson-Newton- 
formula. The right hand side contains Dirac distributions at the frequencies, hence it 
is at least a second derivative of a continuous function. In the left hand side we have 
W(f) that is the d-th derivative of a continuous function. This gives d > 2. 

Also we know that d < o + 1, hence o > 1. □ 

3.3. Symmetric Poisson-Newton formula. Let f(s) be a Dirichlet series with 
exponent of convergence d, and fix o as before. We have defined a distribution 
W(f, o~)(t) = {^2n p e pt ) 1r + supported on R + . If we make the change of variables 
t i— )> —t, we have the distribution W(f,a)(—t) = n p e~ pt ) 1 R _, which is formally 
defined as 

(-l) d ^^((^(-t)-^(0))l K _). 

This is independent of a on R^ and has a contribution at zero dependent on the 
parameter. 

The sum 

W(f,a) = W(f,a)(t) + W(f,a)(-t) 
is a distribution on R, whose only dependence on cr is at zero, and which formally it 
is equal to E P n P eP '*'> * e R - 

Theorem 3.9. For a Dirichlet series f , we have in R, 

9-1 

V?(/, ( x)(t) = 2^c 2 ^f ) + ^ (A,|k|)6| k | 5 <Aik> . 

«=0 keAu(-A) 



W(f, a)(t) = £<A, k)6 k 6 {x , k) 6 S ] » 



Proof. The Poisson-Newton formula for / is 

9-1 

^,^/"k"(A,k) -r / ^ <■ 
keA Z=0 

Making the change of variables t h-> — t, we have 

9-1 

W(/,<r)(-f) = £ < A ' l k DVl <W> + D" 1 )'*^ ' 

ke-A 1=0 
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where |k| = — k, for k G —A. 

Adding the two formulas, we get the symmetric formula stated in the theorem for 

W(f,a)(t) = W(f,a)(t) + W(f,a)(-t) = J> p e^ . 

p 

□ 

Consider a Dirichlet series /(s) = 1 + ^ a n e XnS and let 

f(s) =J^) = 1 + J2 «ne A " s 

be its conjugate. Then / is a Dirichlet series whose zeros are the {p} and rip = n p . 
Also &k(/) = bk(f). The Poisson- Newton formula for / is 

H/(/»(t) = £<A, k)M<A,k» + 

Corollary 3.10. For a rea/ analytic Dirichlet series f , that is f(s) = f(s), we have 
that for a G M., the numbers ci and 6 k are real. 

The converse also holds. 

The last point is due to the fact that the association / \-t W(f) is one-to-one, as 
its inverse is the Laplace transform. 

3.4. Poisson-Newton formula with parameters. Observing that the space of 
Dirichlet series is invariant by the change of variables s h-> as + (3, with a > and 
/3 G C, we get a parameter version of the main theorem. 

Corollary 3.11. Let a > and (3 G C. As distributions on M*^ we have 
e-£« W(/)(t/«)k = ^«(A,k)e-< A ^6 k 5 Q<A , k> . 

keA 

Proof. This results by applying Corollary 13. 6 1 to g(s) = f(as + (3), which is a Dirichlet 
series for a > 0, 

o( s ) = 1 + a n e~ Xnl3 e- (Ka)s . 

n>\ 

The zeros of g are the numbers , and 6 k is changed to e~^ A ' k ^6 k . 

□ 



We can also give a parameter version of the full Poisson-Newton formula. 
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Corollary 3.12. Let a > and (3 e C. ^4s distributions in R we /iawe 

s-i 

e-i* W(/)(f/«) = E c *K)^ fc) + E a < A > k ) e_<A ' k>/3 ^ <W> , 

fc=0 keA 

where the Ck(c[) are the coefficients of the discrepancy polynomial for a — a[ — CT1 ~ 3 ' ?/3 . 
Proof. Just observe that Oi becomes a[ = cri ~ jR/3 for g(s) = /(as + /?). □ 

3.5. Symmetric Poisson-Newton formula with parameters. The space of real 
analytic Dirichlet series is invariant by the change of variables s i— > as + (3, with 
a > and /? 6l, then we get a parameter version of the symmetric Poisson-Newton 
formula of the previous section. 

Theorem 3.13. Let a > and (3 G R. For a rea/ analytic Dirichlet series f , that is 
f{ s ) — f(s) = f(s), we have in R 

e-^mm) + w(f)(-t)) = E^e (p ^ )|t|/a = 

p 

9-1 

= 2E^K)^ 2 °+ E a(A,|k|> e -<W6| k |5 a<A>k> , 

«=0 keAu(-A) 

with a[ = ^. 

In particular, for a — 1 and g — 1 that we use in the applications, we get 

Corollary 3.14. Le£ /3 G R. For a rea/ analytic Dirichlet series f , that is f(s) = 
f(s) = f(s), of genus g = 1 we have in R 

EV^ 1 * 1 = 2^(0-! - fi) 5 + E < A ' |k|)e- (A ' |k|)/3 6| k | 5 <A , k> . 

p keAu(-A) 

3.6. General Poisson-Newton formula. We can observe that in the proof of the 
Poisson-Newton formula we barely used the Dirichlet character of /. The construction 
of G, the discrepancy polynomial Pf, and the Hadamard interpolation /# can be 
carried out in general for any meromorphic function of finite order with its divisor 
contained in a left half plane. The existence of the Newton-Cramer distribution W(f) 
also holds in this generality. The choice of <j\ is only relevant for the structure at 
of the Newton-Cramer distribution and for determining the exponential decay of 
the test functions to which we can apply the distribution formula. Only in order to 
compute the logarithmic derivative /'/ / we have used the Dirichlet series expansion. 
If we take over the proof for a meromorphic funciton / of finite order and with divisor 
contained in a left half plane, we get the following general theorem: 
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Theorem 3.15. Let f be a meromorphic function of finite order with convergence 
exponent d and its divisor contained in a left half plane. Let W(f) be its Newton- 
Cramer distribution and Pf(s) = c +cis+. . .+Cg-\S g ~ l be the discrepancy polynomial. 
We have in R ; 

^(/) = £c^ ) + £- 1 (/'//), 
3=0 

or 

C{WU)) = (W(f), e-°% + = y|| + Pj(s) . 

The inverse Laplace transform ia a well defined distribution of finite order 

when F has polynomial growth on a half plane, which is the case of F — f'/f when / 
is of finite order and has divisor contained on a left half plane. It is defined as follows. 
Take m which is two units more than the order of growth of F, and define 

Ht)= ff^? 1 ^-- 

J -00 (c + iu) m 2tt 
This is a continuous function, which vanishes on R_. It is independent of the choice 
of c (subject to 3fc > <7i). Then 

■= ^m, 

which is a distribution of order at most m — 2. 

More explicitly, for an appropriate test function ip, letting ip(t) = (p(t)e ct , we have 

(£- 1 (F), v .) = (L(t),(-l)V m) W) 

+ °° F(c + iu \-ir<pW(t)e^»^dt 




00 

00 



(C + iu) m ' ' W 27T 



1)m F( C + »«) / /■ W(()e( ^ A* 

(c + iu) m VA / 2vr 

r+°° F( c + iu) r [ {c + iur {t)e(c+lu)tdt \ du 

= J F(c + tu)^(-u) — : 



—00 
+00 



doing m integrations by parts in the penultimate line. 

We can give a symmetric version of the general Poisson- Newton formula: We denote 

W(f,a)=W(f,o-)(t) + W(f,o-)(-t), 
which is a distribution on R. 
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Theorem 3.16. We have in K. 

3-1 

W(f,a)(t) = 2^6^ + {C-\f/m)+C-\flf){-t)). 

1=0 

3.7. General symmetric Poisson-Newton formula with parameters. Let a > 

0, /3 e C. Take a[ = We denote, as a slight abuse of notation, 

e~^W(f, a)(t) = e~^W(f, a)(t) + e^W(f, a)(-t) . 
Theorem 3.17. We have in R 

9-1 

1=0 

Corollary 3.18. For a real analytic function f and j3 G K to the right of all zeroes 
of f , we have in K. 

9-1 

e-^W(f,a)(t) = 2j2c 2l (a[)Si 21 } + e~&C? (2R(/'//)) (*/«)■ 

We prove this using the following result 

Lemma 3.19. For a real analytic function F , and 7 to the right of the zeroes, we 
have 

e-^C-\F)(t) + e^£- 1 (F)(-t) = e^ 1 C: 1 (2MF) (t) 



Proof We have 

/+00 
F(c + iu)e 
-00 

Analogously, 



e ct C'\F){-t)= I F(c + iu)e- tut — 

00 27T 

F(c + iv)e 

2ir 



/+00 
F(c - ivy 
-00 

/_ 



—00 
+00 



Adding both, we get 

. dv 



/+00 
2(&F(c + iu)) 
-00 

as required. □ 



27T 
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3.8. Converse theorem. Conversely we show that any Poisson-Newton formula is 
associated to a Dirichlet series. 

Proposition 3.20. Let Y2 P n pP be a divisor with convergence exponent equal to d 
and contained in a left half plane 9ftp < <j\. Suppose that a Poisson-Newton formula 

, d-2 

E n ^ pt v = E 6 «<v + E c A a) 

' n>l j=0 

holds in WL, in the sense of an equality of distributions, acting on functions with fast 
enough exponential decay. We assume < ii\ < /x 2 < • • ., where (fi n ) are finitely 
many or \i n — > +oo, and J2 n >i \bn\e~ flnCr < oo. for some a < oo. Then there is 
a Dirichlet series f meromorphic on C, whose Poisson-Newton formula is the one 
given. 

Proof. The condition \ n p\ \p\~ d < 00 allows us to define the function G(s,a) as- 
sociated to the divisor, and the corresponding Hadamard interpolation /#(s, er) = 
exp (J G(s,a) ds). The function fn(s,a) is meromorphic on C and has divisor of 
zeros and poles equal to J2 n pP- For the function we have defined a distribution 
W(fn) supported in R + , which equals with J2 n p e ~ pt m ^+ (the precise meaning of 
the later is our definition of the distributional sense given in Section [2]). Therefore, 

W{f H )= (j2n p e-A 1 R+ + . 
V p / j=o 

By our hypothesis, we have the equality 

d-2 

^(^) = E^ + E<*' ) 

n>l j=0 

as distributions paired against functions with fast enough exponential decay. So for 
3?s > a (enlarging a if necessary), we have 

!k = G(a) = (W(f H ),e- ts ) 

JH 

E^ + E4*V ts 

\n>l j=0 

= F(s) + P(s), 
where P(s) is a polynomial of degree at most d — 2, and 



F(s) = /j2b n 6, n ,e- ts )=J2b n e- 



-PnS 

,n>l / n>l 
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Note that the assumptions mean that this can be paired against e~ st for 9fts > a and 
the last sum is uniformly convergent. So for 3?s > a, we have 

f(s):=ef F ^=e-f p ^f H (s). 

Then f(s) is a meromorphic function on C, and in a right half-plane, we have 



f(s) = exp (^6„e-^) = 1 + J2 a n 



n>l 

which is a Dirichlet series. □ 

4. Parameter dependence 

We analyze the dependence of W(f, d', a) on d! and a, which is concentrated at 0. 
We write 

d' 

R A (a,d') = J2^^ A (0,a)5t 1] 
i=i 

for the remainder in §5$). Note that this is not convergent, as the first summand in 
the right hand side §5$) is also not convergent. But 

R A {a Xi d[) - R A (a 2 , d' 2 ) = W A (f, a u d[) - W A (f, a 2 , d' 2 ) 

is convergent. The dependence on d' > d is easy to describe, since the K i A (0, a) do 
converge for I > d. 

Lemma 4.1. Let d' = d + a, a > 1. Then 

a 

W(f, a, d') = W(f, a,d) + J2 ^ 

= W(f, *,d) + d f;( £ ( d ~ I + j ) (-ar^- k K d+j (0, a)) 

k=0 \max(l,k-d+l)<j<a ^ ' J 

In 'particular W(f, a, d) is obtained from W(f, a, d') by removing the higher order 
derivatives of 5q strictly larger than d — 1 and a convergent part in lower order deriva- 
tives. 

Proof. Just note that for n > and a G C, 



k=0 







□ 



22 V. MUNOZ AND R. PEREZ MARCO 

So from now on, we shall restrict to d! — d and write W(f,a) = W(f,a,d). We 
write 

d-l 



R A {a) = W A (f,a) - W A {f,a{) = £(tf I+1)A (0, a)e^ - K l+1A (0,a 1 )e^ 

1=0 

d-i 

= ^r feij4 (cr)^ fc) 



k=0 



with for k — 0, 1, . . . , d — 1 



c M w = E(!.)(-^) , - fc ^ + M(o^) 

z=fc \ ' 
nfe,A(c) = c kjA (cr) - c fe)j4 ((7i) , 

and 

d— 1 

J R(a) = hm J R A ( ( 7) = ^r fc ( ( 7)4 fc) , 

with 



k=0 



r k = limr fciA . 

A 

Note that the last coefficient for k = d — 1 is clearly convergent and the following 
holds: 

Proposition 4.2. We have that r d _i(<r) = iQ(0, <r) - iQ(0, a x ) = E P (^pr - 
iV^(0, <7i) zs a meromorphic function on C wift po/es o/ order d at the p's. 

For k < d — 2, we do not have such an explicit description, but still a similar result 
holds. 

Proposition 4.3. For k = 0,1,..., d — 2, the coefficient rk(cr) is a meromorphic 
function on a e C with poles of order d at the p 's. 

Proof. This results from the parameter dependence of the Hadamard interpolation. 
We have defined earlier 

G(s,a) = -^nJ—--^ =^n 



Then we write 



p~ s (p - a ) l+1 J A ^'" p (p- (r ) dls -p 

/h(s, a) = exp ^ G(s,a) ds^j , 
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and 

Ms, a) f{s) f(s) 

and 

P f (s, a) = c {a) + d(a)s + ... + c d _ 2 (oV~ 2 . 
Assume that is not in the support of the divisor. Then observe that 

taking the derivative with respect to s, j — 0, 1, . . . , d — 2. As G(s, a) has poles of 
order d — 1 at p (as a function on a), the same happens for Cj(a). 

□ 

4.1. Newton-Cramer functions associated to divisors. Consider a divisor D = 
Y2 n pP supported in a left half-plane and of convergence exponent d. This means that 
9ftp < 0"i for some Uj el and J2 PJ to \ n p\ \p\~ d < 00 • Then, for any a G C — {p} and 
d' > d, we have an associated distribution 

W(D,a, d') 

constructed as in Section [2j All these distributions are supported in E + , and they 
are independent of a and d' in K+. Moreover, we denote W(D, a) = W(D,a,d). 

The distributions W(D, a) are Laplace transformable, i.e., they can be paired 
against e~ st , for 3FJs > cx , for some <jq depending on D and a. Define 

g D (s,a) = (W(D,a),e- s % 

and 



f D (s, a) = exp f- J g D (s, a) ds^j . 



Proposition 4.4. All gi)(s,a), for different o, are equal up to the addition of a 
polynomial of degree d — 2 in s. They have simple poles exactly at the p, with residues 



n p . 



The functions frj{s,a) are meromorphic on C. They have the same divisor and 
they differ by a Weierstrass exponential factor of order at most d — 1. 
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5. Applications. 



5.1. Classical Poisson formula. Consider for A = Ai > the entire function 

f(s) = l-e- Xs . 

This function has order 1, exponent of convergence d = 2, genus g — 1, and its zeros 
are 

pk = l ^ k i 



for fceZ. Also we have 



-log/( S ) = -log(l-e- As ) =J2l e 



k=i 

i.e., with our notation A = W, b k = l/k. 

Therefore the Poisson-Newton formula in M*. gives 



oo 

Xks 



\kez / | R . k>i 



The distribution on the left side, when considered without restricting to R+, is even. 
It follows then 

vfcez / |r* fcez* 

Now the distribution on the left side, without restricting to W, is A-periodic on R. 
So we get the full classical Poisson formula, identifying what we get at 0, 



If we start with 



/(*) = 1 - ae~ As , 

for some a 6 C*, we get also the classical Poisson formula. We have for k G Z, 

1, .2tt, 
Pk = jloga + i—k, 

and 

1 

5fc 



k 

Therefore 



(a*/* £ = A E = Aa^ £ 5 Afc , 

\ keZ J |d» fc>l fc>l 
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which gives the same formula as before. 

Application of the symmetric Poisson-Newton formula. As expected from 
the symmetric form of the classical Poisson formula, one can recover the formula in 
a more direct way from the symmetric Poisson-Newton formula. 

It is also interesting to clarify the structure of the Newton-Cramer distribution at 
0. It helps to understand why the Dirac 8 appearing in the right side of the classical 
Poisson formula is of a different nature than the other S\k for k 7^ 0, something that 
was intuitively suspected from the analogy with trace formulas (see a comment on 
this in [Tj, p.2). 

In order to use the symmetric Poisson-Newton formula we compute the discrepancy 
polynomial Pf for f(s) = l — e~ Xs . We have that o\ = is a zero. From the classical 
Hadamard factorization 

sinh(7rs) = us j^J ^1 — —J , 

fcez* 

we get the Hadamard factorization for /, 

f(s) = 2e- As/2 sinh(As/2) = s\e~ Xs/2 JJ ( 1 - — J et . 

Note that this is equivalent to 



GW = i _ E (_!_ _ 1) 



fcez 
A/2 



tanh(As/2) 

Thus Qf(s) = (log A + 27iin) — |s, with n G Z, and 

P f (s) = -Q' f (s) = c = ^ . 

Therefore we can apply the symmetric Poisson-Newton formula (Theorem I3.9P and 
we get 



\s„ + a Y, s 



A A: 

fcGZ* 



k&Z 
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We finally observe that 




fcez fcez 

because we can reorder freely a converging (in the distribution sense) infinite series 
of distributions 



5.2. Newton formulas. We show in this section how the Poisson-Newton formula 
is a generalization to Dirichlet series of Newton formulas which express Newton sums 
of roots of a polynomial equation in terms of its coefficients (or elementary symmetric 
functions). 

Let P(z) = z n + a\z n ~ l + . . . + a n be a polynomial of degree n > 1 with zeros 
«!,...,«„ repeated according to their multiplicity. For each integer m > 1, the 
Newton sums of the roots are the symmetric functions 



From the fundamental theorem on symmetric functions, these Newton sums can be 
expressed polynomially with integer coefficients in terms of elementary symmetric 
functions, i.e., in terms of the coefficients of P. These are the Newton formulas. For 
instance, if for k > 1 






n 





l<il<...<jj.<n 



then we have 



S 



i 



S 2 

s 4 




2S 2 

3S2S1 + 3X3 

4S 2 S? + 4E 3 Ei + 2E^ - 4E 4 
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We recover them applying the Poisson-Newton formula to the finite Dirichlet series 
f( s ) = e - Xns P(e Xs ) = 1 + aie" As + . . . + a n e~ Xns . 
The zeros of / are the (pj y k) with j — 1, . . . , n, k 6 Z, and 

= afe^ k . 

Thus, using the classical Poisson formula, its Newton-Cramer distribution can be 
computed in R as 



n 



p 3=1 k&L 



It 



iiiA 



Now, using the Poisson-Newton formula in ]R+ 

^e^ = ^(A,k)6 k 5 (Aik) , 

p keA 

taking into account the repetitions in the right side, and that A = (Ai, . . . , A n ) = 
(A, 2 A, . . . , nX), we have using the formula ((Tj) for the 6 k 

which gives the explicit Newton relations. Moreover, Newton relations are equivalent 
to the Poisson-Newton formula in W + in this case. 

For example, for m — 4, 

£4 = 4 (6(4,0,0,0) + 6(2,1,0,0) + 6(1,0,1,0) + 6(0,2,0,0) + 6(0,0,0,1)) , 
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and from 

\ 
b, 



(4,0,0,0) 


= is« 

4 1 


'(2,1,0,0) 


= -E?£ 


'(1,0,1,0) 


= S1E3 


'(0,2,0,0) 


2 2 


'(0,0,0,1) 


= -S 4 



we get 



S4 = - 4S 2 S 2 + 4S 3 Sx + 2E 2 , - 4S 4 



5.3. Abel-Plana summation formula. The full Poisson-Newton formula for f(s) 
1 — e~ As is a "half" classical Poisson formula and can be written for A = 1 as 



(9) 



n>0 VfcGZ / 



We check now that this is essentially the Abel-Plana summation formula (see [T], 
and |6j). The Abel-Plana summation formula compares an infinite sum with the 
corresponding integral. 

Theorem 5.1. (Abel-Plana summation formula) Let f be a holomorphic func- 
tion in a domain containing the right half plane H + = {3lz > 0} and continuous in 
H + with 

lim \f(x±iy)\e~ 27rv = 

?/— s-+oo 

uniformly on compact sets of x, and such that 

\f{x + iy)-f{x-iy)\e- 2 ^dy 







exists for x G M+ and tends to when x — > +00. 
We have 

^ (£/(»>- f / w jf /(< ^y * 

and when the sum converges, 

«=o Jo 1 Jo e - 1 
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Proof. We apply the full Poisson-Newton formula (J9j) to a smooth truncation /jy of 
/ with compact support, so in particular /jy is in the Schwarz class. We consider 
/at(£ + iy) = XN(x)f(x + iy) with Xn( x ) equal to 1 in a neighborhood of [0,N] 
and vanishing outside [— e, iV + 1/2 + e] with < e < 1/2. We also require that 
||/1[o,jv] — /atIIl^r) - > when A" — > +oo. The function f N coincides with / and is 
holomorphic in a neighborhood of the strip < 3lz < N + 1/2. Applying to Jn the 
half Poisson formula, we have 

+oo N 



n=0 n=0 



n=0 

-i r+oo 

of(P) + E / fN(t)e 2mkt dt 

i r+oo +°° /"+oo +°° r—oo 

o/(0) + / fN(t)dt + J2 f N (t)e 2mkt dt - E / f N (-t)e 2mkt dt 
Jo k =i Jo k=i Jo 



-N+l/2 +ra / /-AT r—N 

f(t)e 27Tikt dt - / f(-t)e 2nikt dt ) + o(l) 



I pN + 1/2 -roo / «/ 







where o(l) — >■ when A" — > +oo. 

For > 0, in the domain of holomorphy of /, and using Cauchy theorem, we 
decompose each integral into three line integrals over a rectangular contour: 



N riR rN+iR rN 

f(t)e 2nikt dt = / f(z)e 2mkz dz + / f(z)e 2lTikz dz + / f(z)e 2nikz dz 

JO JiR J N+iR 

-N riR r-N+iR r-N 

f(-t)e 2 * ikt dt= / f(-z)e 2nikz dz + / f(-z)e 2mkz dz + / f(-z)e 27Tikz dz 

io ViJJ J-N+iR 

When R — > +oo the second integral of each line tends to because of the first 
hypothesis on /. The substraction of the third integrals gives 

JV r—N rO rO 

f(z)e 2nikz dz- / f(-z)e 27Tikz dz = i / f(N + iu)e~ 27Tku du — i / f(N-tu)e- 2nku du 

N+iR, J-N+iR Jr Jr 

rR 

= -i (f(N + iu) - f(N - iu))e- 27Tku du . 
Jo 

And the second hypothesis on / shows that this last expression tends to when 
A" — > +oo. So in the limit we are left with 
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£ f(n) - [ N+1/2 f{t) dt = 1/(0) + i [ + °° m ~ fi ; U) dt + o(l) , 
n=o Jo 1 Jo e - I 

and the result follows. □ 

5.4. Euler-MacLaurin formula and generalizations. We have from equation (jSJ) 
Proposition 5.2. For m > d we have 



(i-i) 



DP 

i=i 



And passing to the limit in A we have 

Dt 



m j~\m 



The half Poisson formula can be written as 



+oo 1 



n=0 



We prefer to work with g(s) = 1 - in order to make the simplest choice a = 0. 
Observing that 

W(l - e~ s ) = 1 R+ + W n 



we have 

+°° -r /-r — s\ 

n=0 ^ ' 

Now we apply the formula of Proposition 15.21 with o = and d — 1 to <?(s) = — — 
Then we have for j > d + 1 

^ (2ttz^ 

thus lfj(0, 0) = when j is odd. Recalling that 



f (a) = (-ir^ 
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where the (B n ) are the Bernouilli numbers, we have for j even 

K 2l (0, 0) = 2(-l)' (2tt)- 2 ' £ AT 2 ' = 2(-l)'(27r)- 2I C(2Z) = 
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Thus for m even, or by replacing m by 2m 

Theorem 5.3. (Infinite distributional Euler-MacLaurin formula) We have 

+oo m b r> 2m 

5>= 2 5 ° + lK + + E + ^K 2m (t, 0) . 

n=0 1=1 ^ >' 

This strongly reminds Euler-MacLaurin formula in distributional form. We apply 
this to a test function /. Observe that 



^ l - 1 \f) = -f (2l ' 1) (o), 



and 



/ D 2m \ 

(^W) = <*2™./ (am) > 



,.2irikt 



/* - t _ oo 2-irikt I p —2mkt 



Theorem 5.4. (Infinite Euler-MacLaurin formula) For a test function f in 
the Schwarz class 

ra=0 1/0 Z=l V 7 fc>l 1/0 V 7 

In order to recover the standard finite Euler-Maclaurin formula, we substract the 
infinite formula from itself after a translation by an integer N > 0. 

Theorem 5.5. (Distributional Euler-MacLaurin formula) We have 



E*» = V«+ wEt^ 2 '-''-^'^ E 



e 



n=0 



r(2/)r u ,v ' 1 ^ (2**) 



2m 



[ [0,N] 



A particular case is a finite version of the "half" classical Poisson formula that we 
used in the previous section for deriving Abel-Plana formula. 

Theorem 5.6. (Finite half Poisson formula) We have 

(io) f> = \(s + s N ) + i [m + [J2 e2mkt ) Mm ■ 
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Now, we check that this finite distributional Euler-MacLaurin formula gives the 
usual formula. 

Theorem 5.7. For N,m>0 and f a C 2m function in [0, N] we have 

N r N -, m r> 

£/(»)= / fm+\uiv) + fm + Y.^^f {2k - 1) ( N )-f {2k - 1) ^ 

n=0 ^° k=l ^ '' 

+oo „jv 2nikt I -2mkt 

+ (_i)- y / + * — f 2m \t)dt 

Proof. By density of the Schwarz class in the C 2m topology of C 2m functions in [0, N], 
it is enough to prove the formula for a function / in the Schwarz class. Then it is a 
direct application of the previous finite distributional formula. □ 

With our approach we can generalize Euler-MacLaurin formula by choosing the 
parameter cr^O. 

Theorem 5.8. For a ^ we have: 

+oo 1 m ™ 

n=0 i=cf+l 

fc=0 



where 



and 



£=max(d+l,fc+l) 

-) - E p-^by = (2 "»" ( c ('• s?) + ('• -£)) 

where we use Hurwitz zeta function 



+oo 



CM = E 



£(n + g)« ' 

Remark 5.9. The approach presented to the classical Euler-MacLaurin summation 
formula consists in applying our Poisson-Newton formula to the simplest Dirichlet 
series f(s) = 1 — e~ s . It is clear that this admits an infinite number of generalizations 
applying it to arbitrary Dirichlet series. We will develop these questions elsewhere, 
but we want to notice that in general we get summation formulas over the set of 
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frequencies of well behaved Dirichlet series. In particular, we can obtain summation 
formulas for quasi-peroiodic sequences, etc 

5.5. Ramanujan theory of the constant of a series. Ramanujan developed an 
heuristic theory of summation of series from Euler-MacLaurin formula. He associated 
to some divergent series its "constant" that in his words "It is like the center of gravity 
of a body" (see [23J p. 40 of the first notebook). We refer to chapter 6 of [5] for an 
overview of the theory as exposed in Ramanujan's first notebook. As explained by 
Berndt, Hardy attempted to formalize Ramanujan's theory in [T7] in order to make 
proper sense of Ramanujan's constant for summable series. This does not seem quite 
in the spirit of Ramanujan. Its theory makes sense for sums of the type 

+oo 

£/(»), 

71=1 

where the general term f(n) is given by an analytic function / with global properties. 

Writing in the following form our infinite distributional form of Euler-McLaurin, 
theorem 15.31 



n=l 



we see that it is natural to define the Ramanujan class of tests functions / with global 
properties for which 

For these it is natural to define the Ramanujan distribution RC of infinite order and 
the Ramanujan constant of the series as Ramanujan does by 

1 + °° B 

RC(f) = (RC, f) = (0) - £ (^/^(O) • 

In some sense we can write 

+oo r+oo 
^ f(t) dt , 



iK7(/) = X)/(n)- / 

n=l Jo 



and RC(f) weights the "equilibrium" between the divergence of the series and the 
divergence of the integral. 

In this way we can develop rigorously Ramanujan theory, and this will be done 
elsewhere. 
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5.6. Location of the divisor of a Dirichlet series. From the Poisson formula, 
we get that the distribution W(f) is an atomic distribution on Thus the sum of 
exponentials associated to the zeros cannot be a convergent series for t G M*,. But the 
following lemma gives a simple condition which would imply the analytic convergence 
of the sum. For 9\ < 9 2 , 9 2 — 9i < vr, denote by C(6i, 62) the cone of values of s G C 
with 

0i < Arg s < 6 2 ■ 

Lemma 5.10. If {p} C C{6 1 ,6 2 ), then 

W{t) = Y,n P e pt , 
p 

is a holomorphic function in C(ir/2 — Ox , 37r/2 — 62) ■ 

Proof. For t G C(tt/2 - 1; 3?r/2 - 2 ) we have M(pt) < 0, whence 

\e pt \ < 1 , 

and the series Kg(t) defined in Q is holomorphic in that region, so the results follows. 

□ 

From this we obtain the following straightforward corollary: 

Corollary 5.11. The divisor of any Dirichlet series cannot be contained in a cone 
C(0i, 2 ) for vr/2 <6 1 <9 2 < 3tt/2. 

In particular, once established the functional equation for the Riemann zeta func- 
tion and its unique pole at s = 1, we get that it must have an infinite number of 
non-real zeros in the critical strip < 9fts < 1 (this is not a hard result in any case). 

The above motivates the definition of ^-distribution: 

Definition 5.12. A 9 -distribution is a Newton- Cramer distribution W in R, associ- 
ated to a finite order divisor Yl P n pP su °h that the series 

w{t) = J2 n P ePt 

p 

is absolutely convergent and defines an analytic function in compact sets ofM.^_. 

We also talk about germs of ^-distributions for the germ to the right of of 9- 
distributions. Also we define left-directed divisors: 

Definition 5.13. A divisor D = ^2 p n p p is left-directed if it is contained in a left 
cone C(0i,0 2 ) with vr/2 < Q x < 9 2 < 3vr/2. 
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Thus we proved above that no non-constant Dirichlet series has a left-directed 
divisor and: 

Proposition 5.14. If f is a meromorphic function of finite order with a left-directed 
divisor, then its Newton- Cramer distribution W(f) is a 9 -distribution. 

5.7. Hadamard regularization of ^-distributions. For a ^-distribution, one may 
ask naturally what is the relation of the Hadamard regularization at of the germ of 
analytic function W(t) on ]0, e[ and the distribution W. We show that by the local 
Hadamard regularization we recover the local structure of the germ defined by W at 
0. But the theorem we present below has also a global meaning in R + . It shows that 
the Hadamard formula gives a global integral formula for the full distribution W. 

Theorem 5.15. Let W = W(ai,d) be a 9 -distribution in K. associated to a left- 
directed divisor with convergence exponent d > 2. // o~i G M. is the vertex of a left 
cone containing the divisor, then there are constants C > and e > such that for 
t > 

\W{t)\ < Ct- d e^- e)t . 

The distribution W can be computed by Hadamard regularization at for a test func- 
tion ip such that ip = e ait (p G S is in the Schwartz class, 

(W, ip) = W(t)e-^ (W) - m - 1/(0)t - ... - ^^^ 2 ) dt 

C+oc ( d ~ 2 -i /d-l-2 

= / w(t) Ut)-e-^Y.\ E 

Jo V i=o L \ k=0 

/+oo / d-2 \ 



1=0 / 1=0 




where 



and 



ciM = WitYe-^Ra^i^t) dt , 

n u +00 u 

k=0 k=n+l 

Therefore the distribution W only differs by a contribution at from the classical 
Hadamard regularization H, 

d-2 

W = H + J2ci(°i)$o 



>o 



Note that we can extend the integrals to the whole of R since W = in R*_ 
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Proof. From definition 16.11 the divisor is contained in a left strict cone, so by propo- 
sition 15.101 we have that W is holomorphic in a t-cone containing the positive real 
axis, in particular it is analytic in Moreover e ait W is the d-th derivative of the 
uniformly bounded holomorphic function K d in that cone. Since o~% is not part of the 
divisor, we can consider o\ — e instead of o~% and open slightly the cone in order to 
have a left cone with vertex at o~i — e still containing the divisor. Thus we do have 
the better estimate 

\K d (t)\<Ce~^\ 
because we have in the new cone in the t-plane 

K(0 - ai)t) < -eUt < . 



For to £ we have that the distance d(to) to the boundary of the t-cone is 
> C|to| ; and the estimate of the theorem results from Cauchy estimate, 



\W(t )\ = -e^ 



KAz) dz 



c(to,d(to)) y z ~ t( > 



711 



< C\t \~ d e {ai - e)to . 



Note that in the same way we obtain the control of the derivatives, for I = 0, 1, . . . , d 
and t > 



d}_ 

dt l 



{K d {t) - K d (0)) 



dl+1 T , , ^ 



< Ct 



-l-l-et 



Observing that K d (t) — K d (0) — > for t — > in the cone, we obtain that for t — > 0, 



and I = 1, 2, . . . , d — 1 



-^(K d (t) - K d (0)) = o(r l ) 



Define 



m = m - m - - ... - i.v - 2 



(d-2)\ 



and observe that near 



Therefore, for I = 1, 2, . . . , d — 1 we have 



^(K d (t)-K d (0))^^m=0. 



dt d ' 



Using that for t > 



W(t) = e^^- d (K d (t)-K d (0)), 



POISSON-NEWTON FORMULAS AND DIRICHLET SERIES 37 

we have performing d integrations by parts 
(W,p) = (K d (t)-K d (0),(-l) d — d (e^<p)) 

+00 id 

{-lY{K d {t)-K d {Q))—^)dt 





r+00 

I ( 



/ 

Jo 



hoc 



where we used the estimates for the vanishing at the evaluations at and +00. This 
proves the Hadamard regularization formula stated in the theorem. □ 



6. Functional equations. 

We give in this section a precise definition of the property of "having a functional 
equation" . We know no reference in the classical literature. 

We say that the divisor D\ is contained in the divisor D 2 , and denote this by 

if any zero, resp. pole, of D\ is a zero, resp. pole, of D 2 , and \n p (Di)\ < \n p (D 2 )\ for 
all p e C. 

Definition 6.1. The meromorphic function f has a functional equation if there exists 
a* G K. and a divisor D C Div/ contained in a left cone o~\ + C(8i,8 2 ), with tc/2 < 
Q\ < 9 2 < 3n/2, such that Div / — D is infinite and symmetric with respect to the 
vertical line {3?s = a*}. 

Proposition 6.2. /// has a functional equation and the divisor of f is contained in 
a left half plane then a* G M is unique. 

Proof. Otherwise, if they were two distinct values a*, then Div/ would have an 
infinite subdivisor invariant by a real translation and this contradicts the hypothesis 
that the divisor of / is contained in a left half plane. □ 

Proposition 6.3. If f has a functional equation and the divisor of f is contained in 
a left half plane then Div f — D is contained in a vertical strip. The minimal strip 
< !Rs < cr + }, a + < o~\, with this property is the critical strip and a* = cr ~^" CT+ is 
its center. 
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Proof. Since Div(/) has no zeros nor poles for 9fts > o~i, the divisor of Div/ — D is 
contained in a vertical strip due to the symmetry. The minimal vertical strip has to 



Proposition 6.4. If f has afunctional equation and the divisor of f is contained in 
a left half plane then there is a unique minimal divisor D (i.e., with \n p (D)\ minimal 
for all /) 6 Cj, and a unique decomposition D = Dq + D\, Dq and D\ with disjoint 
supports, with D contained in a left cone a* + C(9i, O2), with n/2 < Q\ < 9 2 < 37r/2, 
and D\ a finite divisor contained in the half plane > cr*} ; such that Div / — D is 
infinite and symmetric with respect to the vertical line {3ts = a*}. 

Proof. We start with D minimal as in the definition, and we define D to be the part 
of D to the left of {9fts = a*} and D\ the remaining part. It is easy to see that Do is 
contained in a left cone with vertex at a*. □ 

Theorem 6.5. If f has a functional equation and the divisor of f is contained in a 
left half plane then there exists a meromorphic function x with Div% = D C Div/, 
such that the function g(s) = xi s )f( s ) satisfies the functional equation 



Moreover, we can write x = Xo • R with Div xo = D — t*D\ and Div R = D\ + t*D\, 
where r is the reflexion along Jfe = a* , and R is a unique rational function up to a 
non-zero multiplicative constant. 

The meromorphic function x (orxo)is uniquely determined up to a factor exp h(s — 
a*) where h is an even entire function. If f has convergence exponent d < +00, then 
we can take x of convergence exponent d, and then x is uniquely determined up to 
a factor expP(s — a*) where P is an even polynomial of degree less than d. In 
particular, when f is of order 1 then x an d Xo are uniquely determined up to a non 
zero multiplicative constant. 

Proof. From Proposition 16.31 we know that a* is uniquely determined as the center 
of the critical strip (which is defined only in terms of the divisor of /). Translating 
everything by a* we can assume that o* = 0. By minimality the divisor of x is 
uniquely determined. Then x is uniquely determined up to a factor exp h(s) where h 
is an entire function. If x(s) = (exp h(s))x(s) gives also a functional equation for /, 
then we have 



be compatible with the functional equation, hence a 



2 



□ 



g(2a* -s)= g{s) 




Therefore 



exp(/i(s) — h(—s)) = 1, 



so for some fceZ, 



h(s) — h(—s) = litik . 
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Specializing for s = we get k — and h is even. 

When / is of convergence exponent d < +oo, and since the divisor of x is contained 
in the divisor of /, then we can take x of convergence exponent at most d. □ 

If / is real analytic, then it is easy to see that x must be real analytic up to the 
Weierstrass factor. We will always choose x to be real analytic. Then g — xf is real 
analytic, and we have a four-fold symmetry and g is symmetric with respect to the 
vertical line 9fts = a*. 

It is of interest to dissociate the contribution to the distribution W(f) of the part 
of the divisor of / that comes from the divisor of xo- This part is the ^-distribution 
W(xo) for which we have Hadamard regularization formula. For a test function tp 
such that ip = e ait <f G S is in the Schwartz class, 

(W( X o), <p) = jf" W( X o)(t)e-^ - m - mt "... - J^§-^ 2 ) dt 



1=0 



where 



where 



ci(<Ti) = jjj Wixomt'e-^R^i^t) dt , 

n u +oo i, 

k=0 k=n+l 

Example. For the Riemann zeta function f(s) = ((s) we have a* = 1/2, er_ = 0, 
a+ = 1, D = -2N*, Di = {1}, and 

X ( s ) = vr -^r( s /2) s ( s -i), 

Xo(s) = vr-/ 2 r( S /2), 
P( s ) = s(s-l). 

Note that 

= x(s)C(s) = n- s / 2 T(s/2)s(s - l)C(s) = 2£(s) . 
(using Riemann's classical notation for £). 



Next, we determine when a finite Dirichlet series satisfies a functional equation. 
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Proposition 6.6. A finite Dirichlet series 

N 

?1=1 

satisfies a functional equation if and only if it is of the form 

[(iV-l)/2] 

f( s ) = e" s aM~ M+ * )s + ce (A '"" )s ) , 

i=0 

where c = 1 if N is even, c = ±1 if N is odd. 

Proof. The Dirichlet series f(s) is of order 1. Suppose that there is some xi s ) of order 
1 with zeros and poles in a left cone such that g(s) = x( s )f( s ) i s symmetric with 
respect to 3?s = a*. By translating, we can assume a* = 0. 

The zeros of f(s) lie in a strip, since e~ x " s f(—s) is also a Dirichlet series. There- 
fore x( s ) has finitely many zeros and poles, and hence x( s ) — §^iJ eAts > f° r some 
polynomials Qi(s),Q 2 (s). The functional equation g(s) = g(—s) reads 

N N 

Qi(s)Q 2 (-s) a n e ( ^ Xn)s = Q 2 (s)Q 1 (-s) ^ a n e^> , 

n=0 n=0 

where we have set a = 1, Ao = 0. 

From this it follows that Qi(s)Q 2 (—s) = cQ 2 (s)Q 1 (—s), c G C*. It follows easily 
that c = ±1. Also 0, Ai, . . . , A^r is a sequence symmetric with respect to /i — Ajv/2. 
So Atv-i = 2/x — Aj and ajv_j = a^. 

If iV even, then Aat/ 2 = /i, c = 1, and 

AT JV/2-1 

^ a n e- x " s = e^ s ^(e (_Al+M)s + e (A ^ )s ) + a^e" 8 . 

n=0 i=0 

If A is odd, then 

N {N-l)/2 

J2a n e~ Ks = e^ ^ ai(e ( - Ai+ " )s + ce^'^ 8 ) , 

ra=0 i=0 

where if c = —1, we have x( s ) — se^. □ 

An example without functional equation. Consider the elementary Dirichlet 
series 

(11) f(s) = 1 + aie- AlS + a 2 e" X2S 

with < Ai < A 2 and aj ^ 0. It is an entire function on C of order 1. 
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If Ai,A2 are rationally dependent, then we may write /(s) = 1 + a\ (e As ) fcl + 

a 2 (e As ) fe2 , for Ai = k\X, A 2 = k 2 \, ki, k 2 > and coprime. We can compute the zeros 
solving the algebraic equation 1 + aiX kl + a 2 X k2 = 0. Therefore, the zeros of f(s) lie 
in at most k 2 vertical lines, and they form k 2 arithmetic sequences of the same purely 
imaginary step. 

If Ai, A2 are rationally independent, then we cannot compute the zeros in general. 
We know that they lie in a half-plane 5fe < 04. Also a 2 1 e X2S f(s) converges to 1 for 
9fts — > —00. So the zeros of f(s) are located in a half-plane 9fts > o 2 . Hence in a strip. 
By Corollary 15.111 there are infinitely many zeros in that strip. 

Now, let (p) be the set of zeros. Then A = {k = (£4, k 2 ) e N 2 | (A*, k 2 ) ^ (0, 0)}, 
and 

and 

S ^n p e pt = y^(AiA;i + \2k 2 )b k 5 Xl k 1 +x 2 k 2 , 

k 

on R*. 

By Proposition 16.61 the Dirichlet series (TTT]) does not have a functional equation 
unless A2 = 2Ai. 



6.1. Functional equation and Hadamard factorization. Let / be a meromor- 
phic function of finite order which has its divisor contained in a left half plane and 
which has a functional equation. In order to simplify we assume that o* is not part 
of the divisor. For g(s) = x(s)f(s) we have 

g(2a* -s)= g(s) , 

and when we express this symmetry in the Hadamard factorization, we get 

Q 9jCT *(2cX* - S) = Qg,a*(s) , 

hence if we write 

k 

the symmetry implies that all odd coefficients are zero — a$ — . . . — 0. 
We observe also that 

Qg — Qx + Qf > 
and for the discrepancy polynomials 



P g = P x + P f 
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Also if the exponent of convergence is d — 2, pairing the Weierstrass factors for 
symmetric zeros p and 2a* — p, gives 

Ed-i ( ■ ( 7 ^ = —(s -p)(s- (2a* - p)) . 

\p-a*J d 1 \(2a* - p) -a* J p - a* K HJK V H " 

Therefore for d = 2, the discrepancy polynomial is constant and we have 

co(x) = P x = P xo = c o(Xo) • 
Note that this also holds if we know that g — 2, in particular when d — 3 and / is a 
Dirichlet series by Corollary 13.71 On the other hand the functional equation implies 
that 

c (g) = , 

therefore we obtain: 

Proposition 6.7. Let f be a meromorphic function of exponent of convergence d = 2, 
or g = 2, which has its divisor contained in a left half plane and has a functional 
equation. We assume that a* is not part of the divisor. 

We have 

co(Xo,o"*) + c (f, a*) = • 

6.2. Gauss formula for the logarithmic derivative of the T-function. From 
our general Poisson-Newton formula we get in one stroke Gauss and related formulas 
for the logarithmic derivative of the T-function. 

The T-function was defined by Euler by the limit 

(n - IV 

T(s) = lh ? r-rrr ( i — r^ s > 

hoo s(s + 1) . . . (s + n — 1) 



ii- 



which shows that f(s) = l/r(s) is an entire function with simple zeros at p n = —n 
for n = 0, 1, 2, . . . Hence the exponent of convergence is d — 2. The divisor of / is 
left-directed and the associated Newton-Cramer distribution is the ^-distribution 

+oo ^ 

^(i/r)(t) = = r -^. 

n=0 

From Euler's definition we can derive directly the Hadamard factorization 

1 _ n ,. .■ 

n=l 



,-s/n 



/W = rM =se7 *n( 1 + D e 

v ' n=l 

where 7 is Euler's constant 

,. 1 , f + °° ( 1 1\ , 
7 = hm > - — log n = e dt 
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(see [29J, p. 246 for the integral expression). This gives the Hadamard interpolation 
for oi = 0. Unfortunately o\ = is a zero of / and the Hadamard regularization 
is divergent at +oo. So we choose G\ = 1 (for example) and derive the associated 
Hadamard factorization 



f( 8 ) = J— = — - = e 7(-i) ff i\ " 1 

m r(s) s-ir(s-i) 11 1 1 

Thus we have for some n e Z, 

<5(s) = 27rm — 7 + 7s , 
and the discrepancy is a constant polynomial 

P f (s) = -Q'(s) = - 7 = c . 
Applying the general Poisson-Newton formula to / we have 

c (WU)) = £ = _ 7 r ''»> 



-n 



s-l 
g-n-1 



e"V ' r(s) 

And by the Hadamard regularization formula applied with <j\ = 1, d = 2 and the test 
function = e~ st , = e^ 1- ^', 

/•+oo 1 

= (W(f), e- st ) = / -<T* (eV s < - l) dt 

/■+oo 1 

Therefore we get 

and plugging the integral expression for Euler's constant we finally prove Gauss for- 
mula for the logarithmic derivative of the T-function |29j, p. 246). 

Theorem 6.8. (C.F. Gauss) 

T'(s) f + °° fe-* e~ st 



r(s) Jo V t 



dt . 



So, according to our interpretation, Gauss integral formula is equivalent to all 
Newton relations for the roots e~ n for all exponents t e R+. This interpretations 
seems new. 

In a similar way we can obtain Binet formula for the logarithm of the T-function 
(see the section on the T-function in [29]), as well as general formulas for higher 
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T-functions, the first example being Barnes T-function. This will be developped in 
future versions of this article. 

6.3. Explicit formulas for Riemann zeros. In this section we apply our Poisson- 
Newton formula to the Riemann zeta function. We obtain a non-classical form of the 
Explicit Formula in analytic number theory. The classical forms can be derived from 
our distributional formula. 

Explicit formulas in analytic number theory go back to the original memoir of 
Riemann [23] on the analytic properties of Riemann zeta function where it is the 
central point of the derivation of Riemann's asymptotic formula for the growth of the 
number of primes. It relates prime numbers with non-trivial zeros of Riemann zeta 
function. Despite the mystery about the precise location of the non-trivial zeros, many 
of such formulas were developped at the end of the XlXth century and the beginning 
of the XXth century (see [20]). Later, general explicit formulas were developed by 
A. P. Guinand [13], J. Delsarte [IT], A. Weil [28] and K. Barner [3], these last ones in 
general distributional form. A classical form of this Explicit Formula is the following 
by K. Barner [1] : 

Theorem 6.9. For an appropriate test function (p with Fourier transform (p analytic 
in a large enough strip, we have 

V^( 7 ) =0{i/2) + <p(-i/2) + ±- [ *(t)<p(t) dt 

p,k>l 

where the ^'s run over the non-trivial zeros of (, the p 'sover prim numbers, and 

(t) = - log 7T + U (jr(V 4 + i*/ 2 ) 

Let's see how one can recover this classical formula from our Poisson- Newton for- 
mula. 

We consider the Riemann zeta function defined for > 1 by 

c( a ) = 5>- = zy* l0,IB • 

n>l n>l 

which is a Dirichlet series with A n = log(n + 1) and o\ = 1 in our notation. It has 
a meromorphic extension to the complex plane s G C with a single simple pole at 
s = 1. It has order o = 1 and convergence exponent d = 2. 
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For > 1 we have the Euler product which gives the relation of the zeta function 
with prime numbers, 

c(8) = Y[(i-p-T\ 

P 

where the product is running over the prime numbers p. Thus 

p, fc>l p, i;>l 

The vector of fundamental frequencies is A = (log 2, log 3, log 5, . . .). We have 6k = 
when k has more than one non-zero entry, and 6k = — \jk for (A, k) = fclogp. 

The Riemann zeta function has a functional equation with a* = 1/2, er_ = and 
er+ = 1. We have, using the notations of section El 

9(s) = g(l - s) 
g( s ) = X(s)((s) 
X (s) = 7r- s / 2 T(s/2)s(s-l), 

Xo (s) = 7r-/ 2 r( s /2), 

R(s) = s(s - 1). 

The Riemann zeta function has a single simple pole at p = 1, and simple real zeros 
at p = —2n, for n = 1,2,..., and non-real zeros in the critical strip cr_ = < 3fts < 
1 = <j + , p = 1/2 + 27- The Riemann Hypothesis conjectures that 7 G K, i.e., that all 
non-real zeros have real part 1/2. These non-real zeros are conjectured to be simple. 
Following the tradition we will repeat them according to their multiplicity, so we may 
skip the multiplicities n p = 1 in our subsequent formulas. 

The Riemann zeta function is real analytic, and we can apply the symmetric 
Poisson-Newton formula (theorem I3.9P and the Poisson-Newton formula with pa- 
rameters (theorem I3.13P . More precisely, in order to get the classical formulas and 
exploit the functional equation, we apply the Poisson-Newton formula with parame- 
ters, Corollary EH for p = a* = 1/2, so a' x = 1 - 1/2 = 1/2 and we get 

Theorem 6.10. 

^^-1/2)1*1 = 2cq(C 1/2) 5q _ (log P )p- fc / 2 (5 klogp + (Ulogp) ■ 

The contribution at 0, Cq((, 1/2) is computed in the Appendix, and we have 
Theorem 6.11. 

c (C,l/2) = -^-|-21og2. 
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We can compute explicitly the contribution of the real divisor to the distribution 
on the left handside: 

Wo (t) = - e m +e ~my e '^\t\ = _ e m +e -*m 1 = _ e i*i/2 +e -§iti_l — 

< J 1 — p AM 



l-e- 2 l*l 2sinh|t| 

n>l 1 1 



Note that -e\ l \l 2 W Q = W( X )(t) + W( X )(-t) = W( X o)(t) + W( X o)(-t) + W(R)(t) + 
W(R)(-t). 

So the associated Poisson-Newton formula on R is 

e n|t| + W (t) = 2c (C, 1/2) S - £ (logp)p- fe / 2 (5 klogp + 5_ feIogp ) 



p,k>l 



2C (C, 1/2) 5 - ^ (<Wn + 5- log™) , 

n>l V ™ 



where A is the Von Mangoldt function, A(jr) = logp and A(n) = if n is not the 
power of a prime number. 

For a test function </? e <S in the Schwartz class, consider its Fourier transform 

<p(x) = [ ip{t)e- ixt dt . 
Jr 



Observe that 



0( 7 ) = I v {t)e-^ 1 dt = [ { V {t)e^ 1 + V {-t)e' l ^ )l ) dt . 
Jr Jr + 

By the real analyticity of ((s), the set of non-trivial zeros is real symmetric, (7) 
(—7), hence 



5>(7)= / M*) + ?(-*))[£■ 



Thus applying now our Poisson-Newton formula to the test function <p we get 
${1) + W [<p] =2c (C, 1/2)^(0)- £(logp)p- fe / 2 (^(Hogp) + ^(-Hogp)), 

7 P,k>l 

where W^ofy?] is the functional 

Wfo,] = / W (t)ip(t) dt . 
Jr 

We compute more precisely this functional. We have 
W = -e-^(W( X )(t) + W( X )(-t)) 

= -e-^ 2 (W( Xo )(t) + W( Xo )(-t) + W(R)(t) + W(R)(-t)) . 
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We assume that (p is holomorphic in a neighborhood of the strip \Qt\ < 1/2, then we 
have by the general symmetric Poisson-Newton formula (or by direct computation) 



(-e-W 2 (W(R)(t) + W(R)(-t)),v) = - / e-W 2 (e^+l)(p(t)dt 

Jr 

= - [ 2cosh(\t\/2)(p(t)dt 
2cosh(t/2)ip(t)dt 
<p(t)e t/2 dt- [ ip(t)e~ tl2 dt 



= -<p(i/2) - <p(-i/2) . 

Now, again using the general symmetric Poisson-Newton formula, more precisely, 
Corollary 13.181 with a = 1 and (3 = 1/2 applied to xo that is real analytic, we have 



-e~W 2 (W(xo)(t) + W( Xo )(-t)) = 2c ( X o, 1/2)^ + Cr x % (m 
And using Proposition 16. 7\ 

(12) c (xo,l/2) + c (C,l/2) = 0, 

thus the Poisson-Newton formula applied to an appropriate test function tp is 



m =0(i/2) +fp(-</2) + (c-} 2 (2n 



Now, we have 



- JZ(logp)p fe/2 (v?(fclogp) + <p(-k\ogp)) 

p,k>l 



i. ir'( s /2) 

log 7T + 



SO 



where 



'1/2 





1 _ 


xo(s; 


1 




'Xo 




.Xo 



2 ° 2 r(s/2) ' 
I 2vr J R 



= - log + K f — (1/4 + it/2) 

Thus we recover the classical form of the Explicit formula stated at the beginning 
of the section. Historically this form is due to Barner that gave a new form of the 
Weil functional. Barner's derivation is based on an integral formula, Barner formula, 
that can be directly derived from our general Poisson-Newton formula. 
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Note that our "explicit formula" appears more concise that the classical formula- 
tion, and even more if we use Corollary 13.141 with /3 = 

Theorem 6.12. We have 

^2n p e pw = 2c (C,0) 5 - ^ ( lo SP) (<W P + 5-fcio gP ) , 

P P,k>l 

and 

c (C,0) = -log(27r) . 

We can compute co(C, 0) from the known Hadamard factorization of Riemann zeta 
function. We have (see [26] p.31): 

e bs -r-r e bs -r-r 

CW = 2( s -i)r( 3 /2 + i) n £ -'^> - .(.-Dr^ n^''^ ' 

where the product is over the non-trivial zeros and 

b = log(27r) - 1 - 7 /2 . 

Now, we have 
thus 

c (C, 0) = -Q c = -b - 1 + c (l/r(s/2), 0) . 
But from the Hadamard factorization of the T-function we have 

1 s 



r(s/2) 2 

(13) co(l/r(s/2),0) 



e : i s l[E 1 (s/(-2n)) 



n>l 

thus 

_7 
2 
and 



(14) c (C,0) = -log(27r) . 

We have the final formula: 
Theorem 6.13. 

n p e p W = -21og(27r) 5 - (logp) (S k i ogp + 5_ fc iog P ) • 

P P,k>l 

Remark 6.14. (Newton relations interpretation) Again, by our general inter- 
pretation, the Explicit formula appears as the "Newton relations" which links the 
non-trivial zeros with the primes, the primes playing a similar role than coefficients 
in Newton formulas. 
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Remark 6.15. (General Explicit Formulas) The derivation given of the classical 
distributional Explicit Formula is general and applies to any Dirichlet series of order 1 
with the required conditions. In this sense the Poisson- Newton formula can be seen as 
the general Explicit Formula associated to a Dirichlet series. The structure at needs 
to be computed in general. But when we have a functional equation, one can apply 
the Poisson-Newton formula with the parameter well chosen so that the structure at 
vanishes from the formula (as we have done in the previous section for the Riemann 
zeta function). The divisor on the left cone gives the general "Weil functional" and 
again, by application of the general Poisson-Newton formula with parameters and 
using Hadamard regularization for this ^-distribution we get a general Barner integral 
formula for the functional. Thus we get a general Explicit Formula with the same 
structure as for the classical one for Riemann zeta function. 

6.4. General Guinand equation. The Newton-Cramer distribution W(f) can be 
naturally be decomposed in the form of an hiperfunction W(f) = W+(f) + W-(f) by 
separating zeros with positive and negative imaginary parts. Both W+(f) and W-(f) 
are analytic functions on cones sharing K + in its boundary. If / is real analytic, then 
its zeros are symmetric with respect to the real axes, giving a relation between W + (f) 
and W-(f). This relation plugged into 

w(f) = w + (f) + wM) , 

givesa functional equation for W + (f) which generalizes Guinand functional equation 
for the Cramer function associated to Riemann zeta-function. 

Therefore this proves that we have general Guinand equations for the generalization 
of the Cramer function for general real analytic Dirichlet series. 

6.5. Selberg Trace formula. It is well known that Selberg trace formula was de- 
veloped by analogy with the Explicit Formulas in analytic number theory and that 
this was the original motivation by Selberg (see [25] , [I]). In this section we explain 
this folklore analogy by showing that Selberg Trace Formula results from the Poisson- 
Newton formula applied to Selberg zeta function. The approach is very similar to 
that of the previous section and we have a unified treatment of both formulas. The 
only relevant difference is that Selberg zeta function is of order 2. 

We consider a compact Riemannian surface X of genus h > 2 with a metric of 
constant negative curvature. Let V be the set of primitive geodesies. The Selberg 
zeta function is defined in the half plane 3?s > 1 by the Euler product 

cx( S )=nri( i - er(p)(s+fc) ) > 

peVk>o 

where r(p) is the length of the geodesic p. 
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We have 



-r{p){s+k)l 

I 

p k>0 1>1 



e -r(p)ls 



I 1 - e T ^ 1 

p,i>i 



e -r(p)l/2 1 e ~r(p)ls 



pri l 2siiih(r(p)//2) 



Thus we compute the coefficients 



/ 2sinh(r(p)//2) ' 

and the frequencies 

(\ (p,i)) = -V = t (p) 1 ■ 

One of the fundamental results of the theory is that (x has a meromorphic extension 
to the complex plane of order 2, exponent of convergence d = 3, thus genus g = 2 by 
Corollary 13 .7\ has a functional equation with a* = 1/2, and its zeros are the following 
(see [27], p.129): 

• Trivial zeros at s — —k with k — 0, 1, 2, . . . with multiplicity 2(/i — 1) (2A: + 1). 

• Non-trivial zeros s — 1/2 ± 27„, n = 0, 1, 2, . . ., where 1/4 + 7^ are the eigen- 
values of the positive Laplacian —Ax on X counted with multiplicity. The 
lowest eigenvalue yields two zeros, s = 1 that is simple, and the trivial zero 
s = with multiplicity 2(h — 1) (we exclude the case of 1/4 as eigenvalue). 

For n < write 7„ = — 7_ n . Therefore the Newton-Cramer distribution decomposes 

as 

W{Cx) = V(Cx) + W Q (Cx) , 
where Wq is the contribution of the trivial zeros and V the contribution of the non- 
trivial ones. We compute on IR* with (3 = 1/2 

W (Cx, l/2)(f) = J2 2 ( h ~ X )( 2n + l)e ( ~ n - 1/2)l " 
= A(h - 1) J2( n + l/2)e- {n+1/2m 

n>0 

v d ( 1 
= -A(h- 1) — 



V.2sinh(|t|/2) 
cosh(t/2) 



(sinh(t/2)) 2 
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And we have 

V((x, l/2)(f) = £ C **M = 2^cos( 7 „t) . 

Now we apply the symmetric Poisson- Newton formula with parameter (Corollary 
EZED with p = 1/2, and we get 

Wo(Cx, 1/2) + t/(Cx, 1/2) = 2c (Cx, 1/2) + ^ |(A, (p, 0}|e H < A ' W)>l/2 ^| o (Aj(p , 0) 

o\^7t I / L n COSh(t/2) ^ r(p) r 

V j (sinh(t/2))^ 2 ^4sinh(r( P )|/|/2)^' 

where we used c (Cx, 1/2) = by Proposition 16.71 

This yields the classical Selberg Trace Formula as stated in [7]: 

Theorem 6.16. (Selberg Trace Formula) We have 

V **= n cosh ^/ 2 ) v t (p) r 

^ ^ 2 1 ^ j (sinh(t/2)P + ^4sinh(r(p)|/|/2) 

We can now manipulate the integral expression for the "Weil functional" a la 
Barner, using the general Poisson-Newton formula as we have done in the previous 
section, etc. These computations will be done elsewhere. 

Remark 6.17. (Gutzwiller Trace formula) The Selberg trace formula is just a par- 
ticular case of the Gutzwiller Trace formula in Quantum Chaos (see [14]). We see that 
in general Gutzwiller Trace Formula, that is the central formula in quantum chaos, 
results from the application of the Poisson-Newton formula to the dynanmical zeta 
function of the Dynamical System when this zeta function has an analytic extension 
to the whole complex plane. Thus non-trivial zeros are related to the quantum energy 
levels and the frequencies to the classical periodic orbits. 



6.6. Lifting formulas. The "lifting formulas" developped in this section are exam- 
ples of Poisson-Newton formulas for Dirichlet series of infinite order. They have a 
transalgebraic meaning that will be developped elsewhere. 

We have normalized our Dirichlet series by ao — 1, but we can carry out the same 
analysis in general for 

f(s) = a + ^a n e- XnS , 

n>l 



with ao 7^ 0. 



52 V. MUNOZ AND R. PEREZ MARCO 

We can write 

f(s) = a 

and we have the associated Poisson-Newton formula in 



+ 



p k a o 



where the first sum is over the zeros (p) of /. But we can also write 

/(s) = (ao-l)+ I l + > a n e~ Ks I =(ao-l) + g(s) . 



i) + (i + £' 

V n>l 



Note that the zeros {77} of <? are the preimages by / of a — 1. Hence we have proved 
Proposition 6.18. W^e /lave m where the first sum is taken with multiplicity 

E ^ = 5> k >4w<w 

r,;/(r,)=a -l k fl 

Observe that when ||k|| = 1, say k = (0, . . . , 0, 1, 0, . . .) with 1 at the j-th place, 
then 

b k = -dj . 

Now adding these Poisson-Newton formulas for ao = 1, 2, 3, . . . we get 
Corollary 6.19. In R* + we have 

+oo / +0O \ 

E E ef* + Y,*n°n8x n \= £ (A,k)6 k C(||k||)o <A , k> . 

m =0 \p;f(p)=m n=l J keA;||k||>2 

Or also 

Corollary 6.20. In R* + - {\ n }, we have 

£ et*= (A,k)6 k (2C(||k||)-l)5 (A , k> . 

pe/- 1 ^) k eA;||k|| even 

7. Appendix 

In this appendix we determine the relation between Qf, a and <5/,o- in particular, 
this the variation of the coefficient co(f, a) from co(/, 0) and we apply this to compute 
co(C,l/2). 
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Let / be of finite order and consider the Hadamard factorization of / (see [2] p. 208) 



f(s) = s no e Q f^l[E m (s/p) n 

p^O 



p 



where m = d — 1 > is minimal for the convergence of the product with 

E m (z) = (1- z)e z+1 ^ 



and Qf is a polynomial uniquely defined up to the addition of an integer multiple of 
2ni. Consider now a G C and the corresponding Hadamard factorization centered at 



f(s) = (s- ay-eOt'-M J] E m (j^J 



pj^cj 

We want to understand the difference between these two factorizations. We take 
logarithmic derivatives to get 



— o J ' ^ (p — o) m s — p {—a) m s 

n _ s m 1 s m 1 



Therefore 



s P s — p cr" s — a 



Q'f.a ~ Q'f =no { a \ I °' + n a 



■ M ~ f "" -'-a) m s ' '"° a m {s - a) 

s m (p - a) m - (s - a) m p m 1 



p p m (p -a) m s- p 



For m — 1 this reduces to 



(is) q}, ct - =^ + ^ + £ 

J ' J cr cr p(p — cr 

For m = 2, it becomes 

„, „, 2a - s s + a >r-^ (-2pa + a 2 )s + pa 2 

Qf, ~ Qf = n °^r- + n*—r + L — J(p^o) 2 — • 

We also can calculate the discrepancy when we change from cr to a' by considering 

Qf.a' - Qf.a = (Qf.a' ~ Q f) ~ {Qf.a ~ Q f) ■ 
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For m = 1, it is of the form As + B, where 

rv n a t ^ (a + a' - 2p) (a' - a) 



~~ (rr' - rrV ( rr' - rr^ 



a' -2a a -2a' v-^ (pa + pa' - 2aa')(a' - a) 

n r 



CT v - °? a k - °? ^ p (p - °)\p - °'? ' 

To compute co(C, 1/2), note that Co((, 1/2) = — co(xo, 1/2) by (fl2l . The value of 

c (Xo,0) = -y- + - 

follows from (fl3|) . The zeros of Xo are the negative integers — n, n > 0, and are simple. 
Hence the formula (TT5]) reads (for cr not a pole of Xo) 



-c (Xo, cr) + c (xo, 0) = h 



-cr 



cr ^ (-n)(-n-a) 

n=l ' 

(j < J \ ri n -I- 

r'fcri 



cr V n n + a 

n=l x 



T(a) +7 ' 

where the last formula follows from the expression for the logarithmic derivative of 
the the T- function, the digamma function ip, 

T'(s 



. , , r (s) i ^ / i i 

= WT = 7 + 7 



r(s) s 2 — f \n n + s 

v ' 71=1 v 

which results from its Hadamard factorization. 

Finally, We get 

Theorem 7.1. We have, for a — Z 

( \ 1 if \ 
co(Xo,o-) = — -ip[(T) . 

In particular, for a = 1/2, we have (see [3] entry 6.3.3 p. 258) 

V>(l/2) = -21og2- 7 . 

Theorem 7.2. We have, 

c (C, 1/2) = -co( X o, 1/2) = - | - 2 log2 . 
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